• 13R^CT/PTOD4 SEP 2001 
PCT/GB00/00741 

09/914843 

Mobile Station Position Determination in a Satellite Mobile Telephone Systety a 
Description 

The present invention relates to position determination for a mobile station in a 
satellite mobile telephone system. 

Various proposals have been made for satellite mobile telephone systems. These 
include the Inmarsat-M system, the IRIDIUM™ system described in EP-A-0365885, 
the ICO™ system described in GB-A-2295296 and the ODYSSEY™ system described 
in EP-A-05 10789. 

It is often desirable that a mobile station in a satellite mobile telephone system know 
its own position. It has been proposed to determine the mobile station position in a 
control portion of the network and transmit it to the mobile station. However, it is 
desirable to minimise the transmission of system information to mobile stations in 
order to minimise the power and spectrum use for the communication of control 
information so that these resources are available for the primary function of the 
system which is speech communication. 

According to the present invention, there is provided a method of determining the 
position of a mobile station in a satellite mobile telephone system, the method 
comprising the steps of receiving propagation information from at least one satellite 
at a mobile station, the propagation information being sufficient to define two 
intersecting curves on the earth's surface, and processing the propagation information 
at the mobile station to identify the location of an intersection of said curves for 
determining the position of the mobile station. 

In one form, the propagation information is received from one satellite and comprises 
propagation time data and doppler shift data for the path between the satellite and the 
mobile station. In a preferred embodiment of this form, the method comprises the 
steps of: 

(a) calculating A as hereinafter defined according to equation (6a); 

(b) calculating B as hereinafter defined according to equation (7a); 



WO 00/52495 




PCT/GB00/00741 



-2- 



(c) calculating the latitude of the mobile station according to equation (8); 

(d) calculating values for two possible longitude of the mobile station 
according to equation (9); and 

(e) selecting one of the longitude values in dependence on the cell on 
which the mobile station is camped being east or west of the satellite's 
ground track. 

In another form, the propagation information comprises propagation time data for 
the paths between the mobile station and the satellites. In one preferred embodiment 
of this form, the method comprises the steps of:- 

(a) determining vector v as hereinafter defined according to equation (15); 

(b) determining the co-ordinates of M as hereinafter defined according to 
equations (17) and (18); 

(c) determining r as hereinafter defined and a pair of co-ordinates for P as 
hereinafter defined according to equations (21) and (19); 

In another preferred method of this form, the method comprises the steps of:- 

(a) solving equation (25) for Z as hereinafter defined; 

(b) determining X as hereinafter defined and Y as hereinafter defined for P 
ac hereinafter defined from/(Z) as hereinafter defined and g(Z) as 
hereinafter defined; 

Preferably, the two satellite methods include the steps of converting the cartesian co- 
ordinates of P into longitude and latitude according to equations (22) and (23), and 
selecting one of the longitude and latitude value pairs produced on the basis of the 
identity of the cell on which the mobile station is camped. More preferably, the step 
of determining T pl as hereinafter defined and T p2 as hereinafter defined using equation 
(1) is included. 

In the foregoing, the equation numbers refer to the equations so numbered in the 
following detailed description. 
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The present invention is also embodied in a mobile station for a satellite mobile 
telephone system including transmitting means, receiving means and a controller 
configured for operation in a method according to the present invention. Such 
configuration will preferably largely comprise programming of the mobile station's 
controller to process the propagation information. 

Conveniently, the propagation information and any other necessary data, e.g. satellite 
nadir locations, will be transmitted to a mobile station in a broadcast control channel 
or a control channel associated with the mobile station. 

An embodiment of the present invention will now be described, by way of example, 
with reference to the accompanying drawings, in which:- 
Figure 1 illustrates a satellite mobile telephone system; 
Figure 2 shows a mobile station; 

Figure 3 is a schematic diagram of a satellite access node of the system of Figure 1; 

Figure 4 is a data flow diagram for the system of Figure 1; 

Figure 5 is a plot of propagation time versus distance from satellite nadir; 

Figure 6 is a plot of doppler shift versus distance along satellite ground track from its 

nadir; 

Figure 7 illustrates the geometries of cells served by a satellite; 

Figures 8 illustrate a geometry used for determining mobile station position; 

Figure 9 is a flowchart illustrating a first method according to the present invention; 

Figure 10 illustrates two intersecting planes used for determining mobile station 

position; 

Figure 1 1 illustrates a geometry used for determining mobile station position; 
Figure 12 is a flowchart illustrating a second method according to the present 
invention; and 

Figure 13 is a flowchart illustrating a third method according to the present invention 
Referring to Figure 1, a satellite mobile telephone system comprises a plurality of 
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satellite access nodes la, lb, lc interconnected by a high capacity digital network 2 
(hereinafter "the backbone network"), a plurality of satellites 3a, 3b, a plurality of a 
mobile stations 4a, 4b, gateway mobile satellite switching centres 5a, 5b, 5c providing 
connections between the satellite access nodes la, lb, lc and other networks 6, a 
network management centre 7, a satellite control centre 8 and a tracking, telemetry 
and control station 9. The network management centre 7, the satellite control centre 
8 and the tracking, telemetry and control station 9 are interconnected by a lower 
capacity digital network 10 which is also connected to the backbone network 2. The 
other networks 6 comprise the public switched telephone network (PSTN), cellular 
telephone networks and the like. 

The satellite control centre 8 and the tracking, telemetry and control station 9 control 
the operation of the satellites 3a, 3b, for instance setting transmit power levels and 
transponder input tuning, as directed by the network management centre 7. 
Telemetry signals from the satellites 3a, 3b are received by the tracking, telemetry and 
control station 9 and processed by the satellite control centre 8 to ensure that the 
satellites 3a, 3b are functioning correctly. 

The satellites 3a, 3b repeatedly broadcast control information in broadcast control 
channels (BCCH). A mobile station 4a repeatedly receives and dzzzdes the 
information in the BCCH broadcast by the satellite 3a, 3b currently serving it. 

During a telephone call, a mobile station 4a, 4b communicates with a satellite 3a, 3b 
via a half duplex channel comprising a downlink channel and an uplink channel. The 
channels comprise TDMA time slots on frequencies allocated on initiation of the call 
or re-allocated during a call. 

The satellites 3a, 3b are in non-geostationary orbits and comprise generally 
conventional satellites, such as the known Hughes HS601 model, and may include 
features as disclosed in GB-A-2288913. Each satellite 3a, 3b is arranged to generate an 
array of beams covering a footprint beneath the satellite, each beam including a 
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number of different frequency channels and time slots. 

Referring to Figure 2, a mobile station 4 is generally similar to the units presently 
available for GSM networks and comprises a codec, a controller 16, a microphone 10, 
a loudspeaker 11, a battery 12, a keypad 13, a radio frequency interface, an antenna 
14, a display 15 and subscriber identification module (SIM) smart card. 

The codec comprises a low bit-rate coder, which generates a speech bit stream at 3.6 
kbits/s, together with a channel coder, which applies error correction codes to the 
speech bit stream to produce an encoded bit stream at 4.8 kbits/s. The low bit-rate 
coder is a linear predictive coder. The channel coder uses Viterbi coding. The codec 
also comprises complementary decoders. 

The SIM includes a processor and a non-volatile memory which stores data 
identifying the subscriber and data for use in encrypted communication. 

Referring to Figure 3, a satellite access node 1 comprises a dish antenna 20 mounted 
for tracking satellites, transmitter and receiver circuits 21 including amplifiers, 
multiplexers, demultiplexers and codecs, a visited mobile satellite switching centre 22 
including a controller 23, a visitor location register 24 and a voice mail box unit 25. 
The mobile satellite switching centre 22 is coupled for communications signals to the 
backbone network 2, to a gateway and to the transmitter and receiver circuits 21. 
The controller 23 is coupled for data signals to the visitor location register 24 and the 
voice mail box unit 25 and may also send and receive data signals via the backbone 
network 2. 

The controller 23 responds to addresses on incoming communications signals, from 
the antenna 20, the gateway and the backbone network 2, by controlling the mobile 
satellite switching centre 22 to output the communications signals on the appropriate 
paths to their destinations, i.e. the antenna 20, the gateway or the backbone network 
2. 
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The visitor location register 25 maintains a record of each of the subscribers registered 
with the satellite access node 1. The voice mail box unit 24 provides storage space for 
voice mail messages for subscribers. 

5 

Referring to Figure 4, a database 30, called the home location register, contains 
records relating to each mobile station 4a, 4b. The record contains the mobile 
station's identity (International Mobile Subscriber Identity or IMSI) and the current 
address of the visitor location register 25a, 25b where the mobile station 4a, 4b is 

=rJ0 registered. The visitor location register 25a, 25b contains a part copy of the home 

location register information (e.g. which services are subscribed to), the status of the 

M mobile station (whether it is "local" or "global) the geographical position of the 

rr| mobile station, the address of home location register where the mobile station is a 
subscriber (to enable billing and other data to be collected at a single point), the 

~ 15 currently active satellite access node with which the mobile station is in 

f ij communication via a satellite, an individual enciphering key and the address of an 

1. 1 

!>i associated voice mail box unit location. 

The home location register 30 may be located in the network management centre 7 
20 (see Figure 1) or may be distributed among the satellite access nodes la, lb, lc (see 

Figure 1) or may be in another network, e.g. when a GSM subscriber roams onto the 
satellite network. 

Referring to Figures 1 to 4, a mobile station 4a may be registered with one of two 
25 distinct statuses; "local" in which the mobile station 4a is permitted to communicate 
only through one local area or part of the satellite system network, and "global", 
which entitles the mobile station 4a to communicate through any part of the satellite 
mobile telephone system. 

30 The mobile station 4a performs an automatic registration process (location update), of 
the kind well known in the art of cellular terrestrial communications, on each 
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occasion when the mobile station 4a is used for an outgoing call, when the mobile 
station 4a is switched on, periodically whilst the mobile station 4a is operating and 
when the mobile station 4a has been unpageable for a predetermined period. As is 
conventional, the location update takes the form of transmitting of a signal 
identifying the mobile station 4a (e.g. by transmitting its telephone number on a 
random access channel). The mobile station 4 decides to camp on a particular cell on 
the basis of measurements of signals on beacon frequencies and the random access 
channel used is the one associated with the cell on which the mobile station 4 decides 
to camp. 

The transmitted signal is picked up by one of the satellites 3a. From the received 
random access channel signal, the satellite access node la serving the satellite 3a which 
sets up a dedicated channel for the mobile station 4a. From the random access 
channel signal, the satellite access node la derives the propagation path time and 
doppler shift and communicates this to the mobile station 4a in an immediate 
assignment message. The propagation time is actually send as the difference between 
actual propagation time and a reference propagation time. 

The first message sent on the assigned channel by the mobile station 4a is the "initial 
message" and when this is received the satellite access node la, calculation of the 
mobile station's position using the propagation time and doppler data is triggered. 
The position of the mobile station 4a is added to the "initial message" which is then 
passed on to the visitor location register 25a where the mobile station 4a is registered. 

At this stage, the visitor location register 25a can be used to confirm that the mobile 
station 4a has a subscription for service in the area for which it is located or 
conversely when first registering with the visitor location register 25a to use 
information provided in the initial message to derive the address of the home location 
register and to contact the home location register during which a part copy of the 
home location register information will be transferred to the visitor location register 
25a and the data in the home location register which points to the visitor location 
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register 25a is updated. The location of the mobile station 4a is thus stored in the 
database of the visitor location register 25a. 

As the mobile station 4a and the cells move relative to one another, the visitor 
location register 25a may determine that the mobile station 4a should be handled via 
another satellite access node lb. Consequently, the visitor location register record of 
the mobile station 4a is amended. 

When a mobile terminated call is made to the mobile station 4a, it is directed to the 
visited mobile satellite switching centre 22a at which the mobile station 4a is 
registered. The visited mobile satellite switching centre 22a causes a paging signal to 
be sent to the cells where the called mobile station 4a is expected to be, i.e. via the 
satellites 3a, 3b serving cells covering the mobile station's location. When the mobile 
station 4a receives the paging signal, it enters into a dialogue with the network to 
obtain a channel for the call. If the mobile station 4a does not respond to the paging 
signal, the network will inform the calling party that the mobile station 4a cannot be 
reached. 

When a mobile originating call is to be made, the user of the mobile station 4a enters 
the number to be dialled and presses a "send" key. The mobile station 4r. rerpcnds by 
sending an access request using a random access channel on a dedicated frequency, 
corresponding to a cell, to the network via one visible satellite 3a, 3b. The satellite 
access node la responds on the basis of the cell from which it receives the access 
request. The selected satellite access node la sends a immediate assignment message 
using a access grant channel. The immediate assignment message includes the 
physical parameters of a stand-alone dedicated control channel assigned by the 
network. The radio link setup is completed using the stand-alone dedicated control 
channel and then the connection of the call is attempted. On connection of the call, 
the mobile station 4a is allocated a traffic channel that is then used for the call. 



In a satellite mobile telephone system, as described above, where the satellites are in 
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medium earth orbit, allowance must be made by the mobile stations for different 
propagation times between different locations and a satellite and for doppler shifts 
due to the motion of a satellite. The communication of this information to a mobile 
station will now be described. 

Referring to Figure 5, the propagation time between a mobile station and a satellite is 
a function of the radial distance between the mobile station and the satellite's nadir 
and of elevation. The propagation time is independent of the direction of the mobile 
station from the nadir and can change at up to S.ljas s" 1 . Contours of constant 
propagation time, called Z-arcs, are defined. In most cases, there is only one Z-arc 
crossing a cell. However, the more peripheral cells may be crossed by more than one 
Z-arc. 

Referring to Figure 6, which shows the doppler shift for a rising satellite versus 
distance of mobile station from the satellite's nadir along the satellite's earth track. 
The doppler shift for the uplink can change at up to 4.17Hz s" 1 and that for the 
downlink can change at up to 4.59Hz s' 1 . 

The doppler shift depends on the motion of the satellite relative to the earth's surface. 
For mobile stations on the projected earth track, and crAy there, the latitude 
dependent earth surface velocity and the ±45° angle variations in the compass 
direction of the ground track exactly offset each other. Accordingly, the doppler 
shift varies only with the distance from the nadir. The maximum doppler shift 
occurs on the projected ground track at the horizon. For mobile stations located off 
of the projected ground track, the doppler shift varies with sub-satellite latitude. This 
variation is less than ±lkHz at 2.2GHz for a cell and is deterministic and predictable. 

The satellite produces 163 beams, each defining one cell. There are 19 beam or cell 
types, as shown in Figure 7. 



A mobile station acquires system time and frequency information from the broadcast 



9 



WO 00/52495 




PCT/GB00/00741 



control channel using a combination of energy-profile on broadcast control channel 
signals, fast Fourier transforms on frequency control channel signals and correlation 
with synchronisation data in the broadcast control channel signals. The broadcast 
control channel comprises a burst every 25 timeslots which is transmitted at a higher 
5 power level than traffic channel bursts. The frequency control channel comprises 
two bursts of unmodulated carrier in every 25 timeslots, the frequency control 
channel bursts immediately following a broadcast control channel burst. 

For initial acquisition, the mobile station's receiver tunes between the nominal 
f JO broadcast control channel frequencies used in the system. The in-phase and 
\ s f quadrature components of received signals are sampled at the symbol frequency and 
H then the instantaneous power is calculated from these samples. The power values are 
m averaged over a moving window having a duration of one burst. A broadcast control 

channel burst is identified when the average power exceeds a threshold value. 

l\l When the threshold value has been exceeded, the identification of a broadcast control 

I -i channel burst is checked. First, using the timing of the provisionally identified 

J f broadcast control channel burst as a reference, a window is opened in the middle of 

the period expected to be occupied by the frequency control channel bursts and the 
20 received signal sampled. A fast Fourier transform of these samples is then calculated. 

If the frequency control bursts are present, a component at the carrier frequency, 

allowing for the effects of doppler shift, should appear in the fast Fourier transform 

result. 

25 If a potential carrier is detected in the fast Fourier transform result, the mobile station 
local oscillator can be adjusted exactly to the correct frequency. The broadcast 
control channel signal is then correlated with a synchronisation word to identify a 
reference point in the broadcast control channel signal. 

30 Once the mobile station has synchronised with the broadcast control channel, it can 
decode the information contained in the broadcast control channel. This information 
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includes the identity of the current cell, the location of the satellite's nadir, the 
satellite transmitting the decoded broadcast control channel, the identity of 
neighbouring cells and the broadcast control channel frequencies for those cells. The 
mobile station then tries to detect and synchronise with the broadcast control 
channels of the neighbouring cells. 

The mobile station compares the strengths of the broadcast control channels detected 
and determines that it is in the cell, served by the strongest broadcast control channel 
signal which, if possible, is not on the satellite's ground track. 

Once the mobile station has selected a broadcast control channel, it decodes the 
broadcast control channel information. The information includes the propagation 
time or times (T 0 ) between the satellite and the or each Z-arc crossing the cell. The 
mobile station will generally not be on a Z-arc and allowance must therefore be made 
for variations in actual propagation time (T p ) corresponding to the mobile station's 
position relative to the closest Z-arc. 

After decoding the broadcast control channel information, the mobile station 
transmits a random access channel burst via the satellite to the satellite access node 
serving the satellite. The random access channel has a dedicated frequency different 
from that used for other channel types. On receiving the random access channel 
burst, the satellite access node calculates the propagation time offset (T p - T 0 ) from T 0 . 
The satellite access node then transmits to the mobile station 2(T p - T 0 ) in the 
immediate assignment message in the AGCH (access grant channel). 

The mobile station obtains doppler shift information in a similar manner, a reference 
doppler shift (f 0 ) is decoded from the broadcast control channel and a differential 
value 2(f d - f 0 ) is transmitted to the mobile station in the immediate assignment 
message in the AGCH. 



The mobile station now has the following parameters T 0 , 2(T p - T 0 ), f 0 , 2(f d - the 
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altitude of the satellite, the longitude and latitude of the satellite's nadir together with 
their respective rates of change which are obtained from the broadcast control 
channel and the velocity of the satellite (programmed into the mobile station during 
manufacture). The mobile station can then solve the equations defining the constant 
propagation time and constant doppler shift contours on which it is located to 
produce two values for the intersections of these contours. The intersections will be 
in cells either side of the satellite's ground track, since the mobile station has selected 
a cell not on the satellite's ground track, and the correct intersection can be identified 
from the mobile station's knowledge of the cell's identity and the beam pattern. 

Referring to Figure 8, T p can be calculated by the mobile station from:- 

T P =^T 0 (1) 
where r = 2(T p - T 0 ). 

T p is related to the satellite's latitude <j) s and longitude A. s and the mobile station's 
latitude <|> T and longitude A^- by:- 

T p = — -y/r 5 2 + Tj - 2r s r T [cos <f> s cos <j> T cos(/L 5 - Xj ) + sin <j> s sin <fi T ] (2) 

The doppler shift f d can be calculated by the mobile station from:- 

/,=§ + /<> O) 
where 6 = 2(f d - Q. 

The doppler shift is related to the satellite's latitude <j> s and longitude X, s and the 
mobile station's latitude <|> T and longitude 7^ by:- 
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ca 

- sin <f> s cos <p T cos(/l 5 - ^ )$ s — cos ^ s sin ^ r cos(/L 5 - A r ){i r 

- sin $ s cos ^ r cos(/l s — )(>i 5 — Xj ) 



_ Zk^L [( C0S < f >s s i n ^ ^ + ( S i n ^ cos ^ )^ ] ( 4 ) 



where r s is the distance from the centre of the earth to the satellite, r T is the distance 
from the centre of the earth to the mobile station, c is the velocity of light and d is 
the distance from the mobile station to the satellite, i.e. cT p . 
r r = r s = 0 
$ T = At = 0 



That is, it is assumed that the altitudes of the mobile station and the satellite are 
constant and the mobile station is stationary. Consequently:- 



/„ = - 



r S r T 



cd 



— sin^ s cos^ r cos(A s — Ay — sin <j> s cos<p T cos(/l s - Xj. )X S 

+ (cos^sin^)^ s 



(5) 



From (2) and (5): 



0 



A = 



d 2 -rl-r' T 
~ 1r s r T 



(6a) 



= cos <f> s cos <f> T cos(/ty - ) + sin <p s sin <f* T (6b) 



and 



5 = 



_ fad 



(la) 



= - sin <f> s cos 4> T cos(/L 5 - A T )$ s - sin <p s cos <j> T cos(/l 5 - X T )X S 

+ (cos <f> s sin $ T )$ s (lb) 

From (6) and (7):- 
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^ r = sin 



-i 



a($ s + X s )sin <f> s + B cos $ s 
j> s + X s sin 2 <f> s 



(8) 



^ is in always in the range -7t/2 to 7i/2 and, consequently, no ambiguity arises 
when determining $r when onlysin -1 0 T is known Furthermore, the motion of the 
satellite is such that ^ + X s sin 2 <f> s * s never 0. 

Once, ^ is known, can be obtained from:- 



cos(/l 5 — A r ) = 



f A - sin ^ sin 0 r ^ 
cos </> s cos ^ r 



(9) 



There are two solutions to (9) for , namely ^ anc j 



A r = A 5 — cos" 



^T 2 = ^5 + COS 



-1 



^4 - sin 0 5 sin 0 r 
^ cos <f> s cos<p T j 

f A - sin 0 5 sin <fi T ^ 
^ cos^ 5 cos^ r 



However, it should be borne in mind that cos <|> s is never equal to 0 and if cos0 r e ^ 
then <|> T = 7i/2 or -7c/2 and this would be enough to determine the position of the 
mobile station. 



Accordingly, referring to Figure 9, the mobile station performs the following steps to 
determine the latitude and longitude of its location:- 

(a) calculate T p according to equation (1) (step si); 

(b) calculate f d according to equation (3) (step s2); 

(c) calculate A according to equation (6a) (step s3); 

(d) calculate B according to equation (7a) (step s4); 
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(e) calculate according to equation (8) (step s5); 



(f) determine if $ T = n/2 or -n/2 (step s6) and if not:- 

(i) calculate values for jj^ according to equation (9) (step s7); and 

(ii) select or X T2 in dependence on the serving cell being east or west 



of the satellite's ground track (step s8); 
else 



(iii) determine sole longitude value (step s9). 

The steps may be performed in any order consistent with the required values being 
available for each step. Steps may also be performed in parallel. 

In a second embodiment, in which signals from main and auxiliary satellites 3a, 3b 
can be received, the mobile station 4a can be considered to be simultaneously in two 
cells. In this case, the mobile station 4a obtains T 0 , 2(T pl - T 0 ) values for the cell 
selected for access. Since the timebase is common at the antenna interfaces of 
different satellites 3 a, 3b, it can measure (T pl -T p2 ) and, using the above information 
and equation (1) the mobile station 4a can calculate T p2 . 

The mobile station is at the intersection of three spheres. The first sphere has its 
centre C at the centre of the earth and a radius R equal to the radius of the earth (The 
earth is assumed to be spherical). The second sphere has its centre S a at the main 
satellite 3a and a radius D a equal to the distance from the main satellite 3a to the 
mobile station 4. The third sphere has its centre S b at the auxiliary satellite 3b and a 
radius D b equal to the distance from the auxiliary satellite 3 b to the mobile station 4. 

The mobile station 4 is notified of the latitudes <|> Sa , <|> sb and the longitudes X^, A^ b of 
the main and auxiliary satellite's nadirs in the broadcast control channels of the two 
cells. The latitude and longitude values are converted into earth centred earth fixed 
cartesian co-ordinates:- 
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X s = (R + H) cos <j> s cos /l 5 
7 5 = (J? + 7/ ) cos sin X s 

a a a 

Z s = (R + H) sin </> s 



X S/> = (R + H)cos0 S{> cos/l s 

= (/? + H)cos0 Sb sin 
Z Sb = (R + H) sin 4> Sb 



(10) 



where H is the altitude of the satellites. 



The first sphere is defined by: 



X 2 + Y 2 +Z 2 =R' 



(11) 



From (10) and (11), the following simultaneous equations can be obtained: 

X 2 +Y 2 + Z 2 = R 2 

-2(XX Sa +VY So +ZZ s J = D 2 a -R 2 -(X 2 Sa + Y 2 o + Z 2 Sg ) (12) 
[- 2{XX h + YY Sb + ZZ h ) = D 2 b -R 2 - (X^ + l£ + Z\ ) 



Let:- 



D 2 -R 2 -(X 2 Sb+ Y 2 b+ Z 2 b ) 
~2 



(13) 



Sub (13) into (12):- 

X 2 +Y 2 +Z 2 =R 2 
XX So+ YY So +ZZ So =C a 
XX s . + 1X +ZZ X . =C k 



(14) 



XX Sa + YY So + ZZ Sa = C a defines a first plane P, and XX ^ + YY S + ZZ^ = C b defines 
a second plane P b (Figure 10). 
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v is a vectoralong the lineD whichis the intersection of the planes P aJ P b < 

(15) 



Z,X, =V, 



Z o X o 



X v Y? — X & Y* — 

*a b b s b s a 3 



To determine D, a point and a vector are needed, A point M on the line D can be 
determined from:- 



X v X + Y, Y + Z< Z = C„ 



X s X + Y s Y + Z s Z = C b 



b -b "A 

aX + /3Y + yZ = 0 



MeP a 

M * p b 06) (Figure 11) 

M e plane (0,S a ,S„) 



a, p, and y must be determined 



a 



is a vector normal to OS a and OS b . 
Y 

Then: n 3 = OS a *OS b 
X, 



Y s a Z s b ~ Z s a Y s b 
Z s a X s b ~ Z s h X s 0 



Thus: - 

a ~ Y s a z s b ~ Z s a Y s i 
P = Z s 0 X s b - z s b x s 

r = x s a Y s b ~ Y s a x s b 



(17) 



Consequently, M = (X M> Y M , ZjJ is the solution of three simultaneous equations with 
three unknowns:- 
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08) 



A parametrical equation of D is: 
P = M+rv 



x = x M+rVl 

Y = Y M +rv 2 
Z = Z M +rv 3 



(19) 



Sub (19) into (11):- 

(X M +rv,) 2 +{Y M +rv 2 y +(Z M +rv 3 ) 2 = R 2 



(20) 



Thus:- 



(Vl < + +v » )r a +2(X M v t +Y M v 2 + Z M v 3 )r + X 2 M +Y 2 + Z 2 M - R 2 =0 (21) 



M 



i II Equation (21) is a quadratic equation in r and can be solved because X M , Y M and Z 
1 1 are known from (18), v 1? v 2 and V3 are known from (15) and R is the radius of the 

earth. Solving (21) gives two values for r (r l and rj from which two locations P (P t 

and Pj) can be obtained from (19). 
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The cartesian co-ordinates in which Pj and P 2 are expressed can then be converted 
into longitude and latitude values thus:- 



X p = tan' 1 



( Y \ 



<l> p = tan 



-1 



(22) 
\ 



+ Y 2 . 

p T I p j 



(23) 



The mobile station can determine which set of latitude and longitude values apply on 
the basis of the identity of the cell, served by the main satellite 3a, on which it is 
camped and knowledge the serving beam. 
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Accordingly, referring to Figure 12, the mobile station performs the following steps 
to determine the latitude and longitude of its location:- 

(a) determine T pl and T p2 using equation (1) (step sll); 

(b) determine vector v according to equation (15) (step sl2); 

(c) determine the co-ordinates of M according to equations (17) and (18) (step si 3); 

(d) determine r and P according to equations (21) and (19) (step sl4); 

(e) convert cartesian co-ordinates of P into longitude and latitude according to 
equations (22) and (23) (step sl5); and 

(f) select one of the longitude and latitude value pairs produced in step si 5 on the 
basis of the identity of the cell on which the mobile station is camped (step 
sl6). 

The steps may be performed in any order consistent with the required values being 
available for each step. Steps may also be performed in parallel. 

In a third embodiment, in which signals from main and auxiliary satellites 3a, 3b can 
be received, the mobile station 4a can be considered to be simultaneously in two cells. 
In this case, the mobile station 4a obtains T 0 , 2(T p - T 0 ) values for both cells and then 
determines the intersections of the contours of constant propagation delay for both 
satellites 3a, 3b. 

The values for T pl and T p2 for the two satellites can be obtained from equation (1). 

The mobile station is at the intersection of three spheres. The first sphere has its 
centre C at the centre of the earth and a radius R equal to the radius of the earth (The 
earth is assumed to be spherical). The second sphere has its centre S a at the main 
satellite 3a and a radius D a equal to the distance from the main satellite 3a to the 
mobile station 4. The third sphere has its centre S b at the auxiliary satellite 3b and a 
radius D b equal to the distance from the auxiliary satellite 3b to the mobile station 4. 

The mobile station 4 is notified of the latitudes <|> Sa , <f> Sb and the longitudes X $2 , k Sb of 
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the main and auxiliary satellites in the broadcast control channels of the two cells. 



Starting from equations (14), the second and third equations can be used to derive 
expressions for X and Y as functions of Z (i.e. X=/(Z) and Y=g(Z)), on the 
assumption that the determinate of 



*b s b 



*0 



(24) 



(If this determinant equals zero the Y and Z can be expressed as functions of X or X 
and Z can be expressed as functions of Y.) 



C a Y s -CJ S -Z(Z S Y s -Z S Y S ) 

X = f(Z) = 

Xs o Ys b ~ X h Ys o 

Y = g (2) = CbXs ° ~ C " Xs t ~ Z ( Z s„ X s a - Z s a Xs b ) 

~ Ys a Xs b +Y S b Xs o 



Substituting these functions into the first equation at (14):- 

(/(Z)) 2 +(g(Z)) 2 +Z 2 =7? 2 (25) 



Equation (25) is a quadratic in Z and can be solved to produce two values of Z (Zj and 
ZJ. From these two sets of co-ordinates for P: 

P 1 = (X 1 =/(Z 1 ) > Y 1 =g(Z 1 ), ZJ and P^pC.^Z,), Y 2 =g(Z 2 ), Z^ 



The longitude and latitude values for P a and P 2 can then be found using equations (22) 
and (23). 

The mobile station can determine which set of latitude and longitude values apply on 
the basis of the identity of the cell, served by the main satellite 3a, on which it is 
camped and knowledge of whether that cell is north or south of the main satellite's 
ground track. 
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Accordingly, referring to Figure 13, the mobile station performs the following steps 
to determine the latitude and longitude of its location:- 

(a) determine T pl and T p2 using equation (1) (step s21); 

(b) solve equation (25) for Z (step s22); 

(c) determine X and Y for P from/(Z) and g(Z) (step s23); 

(d) convert cartesian co-ordinates of P into longitude and latitude according to 
equations (22) and (23) (step s24); and 

(e) select one of the longitude and latitude value pairs produced in step s24 on the 
basis of the identity of the cell on which the mobile station is camped (step 
s25). 

It will be appreciated that many modifications may be made to the above-described 
embodiment while still achieving the same aim. Since the satellites follow predicable 
orbits, information concerning the position and movement of the nadirs of satellites 
need not be broadcast to mobile station but may instead be programmed into mobile 
station, at manufacturer of the mobile station or the user's SIM. The programming 
may take the form of look up tables or algorithms relating time and satellite identity 
to nadir position. Thus, de minimis, the only information that it is necessary to 
transmit from a satellite to a mobile station is the propagation data peculiar to the 
existing positional relationship between the satellite and the mobile station, 



